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We have shown that a local-density-dependent term, 
F[p{r)] p{r), of the pion-nucleus potential with a nuclear den- 
sity, p(r), can be transposed into a conventional density-linear 
term, F[pe) p{r), with an effective nuclear density, pe, which 
is found to be close to ~ 0.6 p(0) for most 7r~ bound states. 
The presently found relations, generalized Seki-Masutani re- 
lations, for the density-quadratic term, the medium-modified 
isovector term and the double-scattering isoscalar term of the 
s-wave pion-nucleus interaction assure that the constant pa- 
rameters in the conventional pion-nucleus potential are inter- 
preted as being effective ones in the light of density-dependent 
effects. 



1. Introduction 

The so called "anomalous s-wave repulsion" in the 
pion-nucleus interaction has been a long-standing prob- 
lem in pion physics (see, for instance, [1]). Many authors 
(see, for instance, [1,2]) have claimed the existence of a 
large extra repulsion arising from the density quadratic 
term, ReBo p{r)'^, from global fits of pionic atom data to 
the Ericson-Ericson potential [3,4]: 



V{r)=U,{r) + Up{r), 
Us{r) = -{2TT/m^) 

X [ei{bap{r) + biAp{r)) + e2Bop{rf] 
Pir) = Pp{r)+Pn{r), 
Ap{r) = pn{r) - pp{r), 



(1) 



(2) 
(3) 



where Us{r) and Up{r) are for s-wave and p-wave parts, 
respectively, and Pp{r) and Pn{f) are the proton and 
neutron density distributions, respectively, and ei — 
1 -|- rriT^/MM and £2 = 1 + m^/2Mjv (hereafter, we adopt 
units of m~^, and m~'^ for &o> and Rei3o, re- 

spectively). New light has been shed on this problem, 
namely, the possibihty of an anomalous repulsion in the 
isovector part (6i) as a unique signature of chiral symme- 
try restoration was pointed out [5,6]. In this respect, the 
deeply bound 7r~ states which were produced in heavy 
nuclei (in 207pb [7-10] and in ^ospb [11,12]) are of par- 
ticular importance, since the binding energies and widths 
of Is 7r~ states depend nearly entirely on the s-wave po- 
tential {Us) and, thus, only the Is states in heavy nuclei 
[N > Z) provide unique information on the isovector 
part of the s-wave potential. In fact, the above con- 
jecture has been proved by using new information from 



the Is 7r~ state in ^°Pb. The experiment indicates an 
enhanced isovector parameter, which is related to a de- 
creased order parameter of chiral symmetry breaking in 
the nuclear medium, 6("76i = ff/f^:^ = 0.781:°;^^ [12]. 
Recent global fits of pionic atom data by Friedman [13,14] 
have shown consistency with the above view. 

Since the potential parameter hi is nuclear-density de- 
pendent in the theoretical context of Weise [5] , a ques- 
tion can be raised as to the mutual relation in adopt- 
ing the conventional Ericson-Ericson potential with con- 
stant coefhcients (model C) and a modified local-density- 
dependcnt (model LDD) potential. In the present note 
we show that both methods are equally valid and are 
connected to each other, since a density-dependent term 
can be effectively linearized. Our procedure gives an 
answer to the following question: where do the bound 
pions probe the nuclear potential? This problem is 
closely related to the so-called Seki-Masutani (SM) re- 
lation [16], which emphasized a strong correlation be- 
tween the density-linear isoscalar parameter (&o) and the 
density-quadratic parameter (Rei?o)- In a sense, we gen- 
eralize the SM relation. 



2. Where does a bound tt 
tential? 



probe the nuclear po- 



First, we examine which part (r) of the nuclear den- 
sity, pir), is probed by a bound 7r~. Two typical nuclei, 
^^O and ^"^Pb, are considered. The former possesses 
only shallow bound states, whereas the latter accommo- 
dates deeply bound states. The nuclear matter densi- 
ties are assumed to take a 2-parameter Fermi distribu- 
tion with a half-density radius (Cp and c^) and a diffuse- 
ness parameter {zp and z„). The parameter values were 
chosen as follows: Cp = 6.652 fm, c„ = 6.892 fm, and 
flp = a„ ^ 0.5234 fm for ^ospb, and Cp = c„ = 2.482 fm 
and ap = a„ = 0.5234 fm for ^^O. The central nucleon 
densities are: po = 0.150 for ^"Spb and 0.172 for ^^^O. We 
take known potential parameters for the p-wave parts [1] 
and a set of s-wave parameters (6o = —0.028, hi = —0.12, 
ReBo = and ImSo = 0.055). 

Figure 1 shows the 7r~ densities {Rni{r)^), the nuclear 
densities (p(r)) and the overlapping densities (namely, 
the nuclear densities probed by tt^), defined as 



S{r)=p{r)\Rniir)\''r^ 



(4) 



for (Is and 2p) and ^ospb (Is, 2p and 3d). From 
these figures we notice that the overlapping density is 



1 



peaked at a radius slightly less than the half-density ra- 
dius, nearly independent o/the nucleus and the tt^ quan- 
tum numbers. This means that the boimd tt^ effectively 
probes a fraction of the full nuclear density (po = p(0)): 



0.60 po- 



(5) 



This is a key to intuitively understanding the following 
results of numerical calculations. 

3. Seki-Masutani relation 

Seki and Masutani [16] emphasized the presence of a 
strong correlation between h{) and ReBo in describing 
pionic-atom binding energies and widths, which we call 
the Seki-Masutani relation. Toki et al. found the same 
correlation also exists for the deeply bound Is and 2p 
states in ^"^Pb [17,18]. Recently, this correlation has 
been revisited, and the following common relation has 
been established both theoretically and empirically [12]: 



60 + 0.215 ReBo = constant. 



(6) 



This means that the binding energies (and widths) are 
nearly unchanged by varying either bo or ReBo as long as 
these parameters are moved together to fulfill the above 
SM relation. In other words, none of 60 and HcBq can 
be uniquely determined from given binding energies; only 
their combination, 6q, can be experimentally determined. 

Actually, various global fits gave widely distributed 
values with large errors for these parameters. For in- 
stance, the best-fit values of {bo, Rei3o} are {0.024(15), 
-0.26(3)} in KLTK90 [2] and {0.000(15), -0.14(7)} in 
BFG97 [1]. With a different definition of 60, namely, 
^0 = ^0 — A^o^, in which the double-scattering term 
(A6o^, to be discussed later) is subtracted, Friedman 
obtained three different sets for {bo, ReBo}: {0.018(19), 
-0.14(4)} (F02a) [13], {0.030(10), -0.21(4)} (F02b) [14] 
and {0.020(10), -0.15(4)} [14]. These widely split val- 
ues are shown by large open circles with both vertical 
and horizontal error bars in Fig. 2. They appear to be 
almost meaningless, contrary to their claim that the pa- 
rameter ReBo is well determined. On the other hands, a 
fit with a gridding ReBg yields a series of points, as 
shown by closed squares (FG98 [15]) and by small open 
circles (F02a [13]), both of which indicate linear relations 
parallel to the SM correlation. 

Recently, the value of bo with a gridding ReBg was 
precisely determined by fits of the Is tt^ binding en- 
ergies in six light symmetric nuclei (^^C, ^^N, ^^O, ^"Ne, 
2^Mg and ^^Si) [12]. The obtained data sets, shown by 
closed circles in Fig. 2 (their sizes correspond to the error 
bars), prove the above SM relation (6) perfectly, yielding 
a precise value of 



60 = -0.028 ±0.001. 



It is to be noted that the bg composed of the widely 
distributed values of bo and ReBo in each of KLTK90, 



BFG97, F02a and F02b would yield a value, close to 
—0.030 ^ —0.034 (after subtraction of the double scat- 
tering term in the case of F02a and F02b) , which is not 
so much distributed as the individual values of 60 and 
ReBo- These authors seemed to be puzzled by the large- 
ness of their |ReBo| values compared with ImBo = 0.055. 
On the other hand, they never took seriously the other 
discrepancy that their 60 's (after subtraction of the dou- 
ble scattering term) take large positive values, meaning 
a strongly attractive 60 in contrast to the nearly vanish- 
ing free ttN value [19,20]. Such odd consequences from 
the global fits must be artifacts, arising from regarding 
bo and ReBo as independent uncorrelated parameters. 

4. Seki-Masutani Ansatz and its generalization 

Seki and Masutani [16] showed that the correlation 
between 60 and ReBo can be understood by imposing 
an effective replacement of the expectation value of the 
density-quadratic term as, 



< p{rf > ^ pe < p{r) >, 



(8) 



which is expected to hold for any pionic atom state with 
a common value of pe (Seki-Masutani Ansatz). They 
claimed an effective nuclear density pe ~ 0.50 po- We 
understand that the localization of the overlapping den- 
sities, S{r), near r ^ Tc, as shown in Fig. 1, is the key to 
justify the SM Ansatz. 

We now examine how the Seki-Masutani Ansatz 
applies to various functional forms of the potential, 
Us{p{r)). Let us add the following form to the real part 
of the s-wave pion-nucleus potential: 



A[/,(r) = F[p{r)\p{r), 



(9) 



where F[p{r)\ is a local-density-dependent (LDD) func- 
tional coefficient, representing an additional density- non- 
linear potential involving a local density, p(r). We solve 
the Klein-Gordon equation with LDD potentials and 
compare the numerical results with those without invok- 
ing LDD assumptions, thereby examining how the effects 
of the LDD potentials are expressed. Following the spirit 
of SM, we study how this effect is transposed into an 
effective coefficient. 



F[p{r)]p{r)^F{pe)p{r), 



(10) 



where -F'(pe) is an effective constant parameter involving 
an effective nuclear density, pe. Namely, 



npe) = 



< F[p{r)] pjr) > 
< P{r) > 



(11) 



(7) 5. Polynomial case 



Let us first consider the case of a density-polynomial 
term, which can be expressed as 
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ReC/«(r) = -(27r/m^)ei[6 + Bp(r) + Dp{rf]p{r), (12) 

where 6 = 60 + — Z)/A x 61 . We performed numer- 
ical calculations of the Is binding energies of ^^O and 
^°^Pb, and found strong correlations between B and b 
and between D and h, as expressed by 



2°^Pb: 6 + 0.25B + 0.031£) = const, 
^^O : 6 + 0.20B + 0.025£) = const. 



(13) 
(14) 



The first relation is the well known Scki-Masutani rela- 
tion for the Bp{r)'^ term, and the second relation is a 
new relation for the Dp{r)^ term. 

The effect of the density-polynomial potential is as- 
serted to be equivalent to the effect of an effective poten- 
tial. 



UsM{r) = -(27r/m^)ei 
x[b + 2.826Bp('^ 



=2) 



(15) 

8.000D(p("=3))2]p(r), 



with pe in units of fm ^. We obtained the following ef- 
fective densities: 



,("=2) 
o("=3) 



for Pb, and 



^("=2) 

n(n=3) 



0.088 . 
0.062 . 



0.071 
0.056 



0.59po, 
0.41po, 



0.42po, 
0.33po, 



(16) 
(17) 



(18) 
(19) 



for ^^O. The effective density for the n = 2 term is around 
0.5 po, as given by SM, but that for the n = 3 term is 
smaller (the pion probes the p{r)^ component more at a 
larger r). 

The effective density does not depend much on the 
individual states, as we have found in the cases of ^''O 
and 2°8pb. We are therefore tempted to expect that the 
concept of the effective density applies to any functional 
form of U{p{r)). 

6. Density-dependent 61 pareimeter 

Next, let us consider the case of an isovector s-wave po- 
tential which is expected to be enhanced due to a partial 
restoration of the chiral symmetry in the nuclear medium 
[5,6]. We assume the following potential according to 
Weise [5]: 



6i(r)Ap(r) 



b{ 



free 



1 — a p{r) 



Ap(r), 



(20) 



where a is a constant parameter representing the density- 
dependent effect. We examine whether or not the effect 
of the above LDD potential can be replaced by a conven- 
tional potential with a constant parameter involving an 
effective density pe, as in 



6iAp(r) = ^ \ A/9(r). 



1 



ape 



(21) 



Numerical calculations of the binding energies yielded 
the dependence of B\s on the parameter a, as shown 
in Fig. 3. On the other hand, the standard procedure 
assuming b\ to be a constant gave another relation be- 
tween B\s and 61. The two relations are found to be 
nearly identical, if we take 61 = 61 with 



0.090 ~ 0.60 po- 



(22) 



The above finding indicates that the computational re- 
sults starting from the LDD potential invoking a p{r)- 
dependent bi parameter are equivalent to the results with 
a constant-parameter conventional potential; the LDD 
potential can be safely replaced by a density-linear con- 
ventional potential with a coefficient, eq.(21). 

7. Double-scattering term 

The isoscalar s-wave parameter in the free ttN scat- 
tering, fefl®®, was determined experimentally to be nearly 
zero [19,20] in accordance with the Tomozawa- Weinberg 
theorem [21,22]. On the other hand, the double- 
scattering gives a major correction to the isoscalar term, 
as givenby [3,4,23] 



-[(6/-T + 26?]<^> 
-[(6o/'--)2 + 26?]|-fc^(p), 



(23) 
(24) 



where kplp) = [{Stt^ /2)p{r)]^/^ . This correction is 
local-density dependent, since the correlation function, 
< 1/r > , is density-dependent and 61 can also be den- 
sity dependent. 

Let us first examine how the above relation changes 
when we take into account the LDD effect on < 1/r >, 
while keeping bi as a constant parameter. We solved the 
KG equation with Ab^^ = — /3(p(r))^/^ bf as an isoscalar 
part, and obtained the Big values in ^®0 and ^^^p^-, a,s 
functions of (3, which were compared with Big calculated 

— DS 

with A6q as a constant. Equating the two results in the 
form 



aC = -/3(Pe)^/n?, 



(25) 



we find the effective density, pe, is nearly a constant, 
around ~ 0.40 po for ^^O and ~ 0.54 po for ^°^Ph. Thus, 
the coefficient in the double-scattering correction is by 
a factor of (pe/po)^/^ 0.73 and - 0.81 (for ^^O and 
208p^^ respectively) changed over the conventional for- 
mula with pe = Po, yielding 



Ab^ 



-1.52 6f. 



(26) 



Next, we evaluated the full LDD double scattering ef- 
fect with a density-dependent bi{p). The same form as 



3 



eq.(20) with a as a running parameter was adopted. The 
results were compared with those calculated with bi as 
a constant, as shown in Table I. The effective density, 
Pe/po ~ 0.60, for and ^ospb agre_es well with that 
found for bi{p). This means that the bi effectively con- 
verted from the isovcctor LDD potential term has the 
same effect as the density-dependent 61 (p) in the double- 
scattering formula, eq.(26). 

If we assert that the isoscalar potential strength ^q, 
eq.(7), is composed of the free 60 (essentially zero), the 
double scattering term and ReSo, we can deduce ReSo 
empirically as [12] 

b* — Ab^^ 

ReBo = -2 5 0.038 ± 0.025. (27) 

0.215 ^ ^ 

This value looks reasonable, when compared with ImBo 
= 0.055. 

8. RIA correction 

The present procedure can be applied to any kind of 
additional potential term. Friedman [13,14] takes into 
account the following correction 

At/RiA = -^l^Mi^+apf - l]kF{pfp{r), (28) 

where do = —0.19 and a = 2.7 fm'^. The effect 

of adding AJ7ria was calculated and was found to be 
equivalent to invoking an additional isoscalar term 

Ab^^^ - -0.021. (29) 

This correction can be obtained also by setting p{r) — > 

pe ~ 0.5/90- 

9. Angle transformation term 

Finally, we make some comments on the so called angle 
transformation term as an additional correction to the 
EE potential, 

A?7at = -(27r/m„)[(m^/2M)erV[cop + ciAp] (30) 
+ {m„/M)e^^V'^{CoPnPp)], 

which has been taken into account in the global fits of 
BFG97, FG98, F02a and F02b, but not in the present 
paper nor in Ref. [12]. This correction behaves like an 
s-wave potential, though it originates from the p-wave pa- 
rameters. We examined how this correction is expressed 
in terms of the s-wave parameters. We found that the 
inclusion of this correction brings the following changes: 

Ab^'^ « +0.003, (31) 

for both and ^o^Pb, and 

AlmS^"^ « -0.008 for^^O, (32) 
w -0.014 for ^°^Pb. (33) 



Namely, this correction brings an attractive effect to the 
real part and a narrowing effect to the imaginary part. 
Thus, the best-fit parameters (ftg) with and without this 
correction have differences, as given above. A small dif- 
ference in 60 between the above global fits and the present 
analysis is attributed to this correction. On the other 
hand, there is no sizable effect on bi: Abi'^ « -t-0.003. 

10. Conclusion 

We have shown that the Seki-Masutani relation be- 
tween &o and HcBq holds for the Is binding energies of 
light and heavy pionic atoms, both theoretically and em- 
pirically. If we admit this fact, the parameter for the 
isoscalar s-wave interaction is well represented by another 
parameter, bQ, as given in eq.(6,7), which takes care of 
the strong correlation between bo and KcBq . The reason 
for this relation is understood from the fact that the over- 
lapping density of any bound 7t~ with the nuclear density 
is peaked at the radius which is slightly smaller than the 
half-density radius; the nuclear density effectively probed 
by a bound tt" is around 0.6 po- 

We have shown that the Seki-Masutani Ansatz can 
be generalized for any functional form of the density- 
dependent potential: a LDD term, F[p{r)]p{r), can be 
transposed into a corresponding term F[pe] p{r) with a 
constant parameter, where pe is an effective parameter. 
Some important cases are summarized in Table II. In 
other words, we can say that the conventional poten- 
tial form persists to be valid, because it can represent a 
LDD potential; the constant parameters in the conven- 
tional potential are regarded as being density-dependent 
parameters at p ~ pe ~ 0.60 po- This also means that 
the observed enhancement of bi (interpreted as a re- 
duction of the chiral order paramater) in ^°^Pb [12], 
R = bl^'^ybi = ff/U'^ = 0.78t°:J^, is for the effective 
nuclear density of p,, ^ 0.60 po- and thus, the reduction 
that would occur for the full nuclear density is expected 
to be - 0.63. 
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TABLE I. Relation between the coefficient a representing 
the density-dependent 6i(p(r)) = bi^°/[l — ap{r)] and the 
equivalent isovector parameter bi in giving the same Is bind- 
ing energies. 



Nucleus 


a 


1 - ^{'■"V&i 


Pe 


Pe/PO 




2 


0.215 


0.1075 


0.63 




3 


0.323 


0.1077 


0.63 




4 


0.428 


0.1070 


0.63 




2 


0.185 


0.0925 


0.62 




3 


0.277 


0.0923 


0.62 




4 


0.370 


0.0925 


0.62 



TABLE IL Summary of the effective densities, Pe/po, ob- 



tained for various types of local-density-dependent functions, 

F[p{r)]p{r). 



F[P] 


Is i«0 


Is '°«Pb 


Remark 


P 


0.42 


0.59 


SM relation 


1/(1 - ap) 




0.60 


density-dependent bi 




0.40 


0.54 


effective Fermi momentum 


pi/V(l-ap)^ 0.63 


0.62 


double-scattering term 
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FIG. 1. Overlapping probabilities (lower frame) of the tt" densities (upper frame) with the nucleon densities (middle frame) 
in typical pionic bound states: (Left) ^'^O. (Right) ^''^Pb. The vertical broken lines show the half-density proton radii and the 
vertical dash-dotted line is for the half-density neutron radius in ^"^Pb. 
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_^ ^ 2 _^ ^ 05 conventional potential with bi as a constant parameter. 

Re Bq 

FIG. 2. Seki-Masutani relation between feo and KeBo- 
Best-fit values of bo versus KeBo as a gridding variable, ob- 
tained in minimization using the Is pionic atom data in 
6 N = Z nuclei, are shown by closed circles, whose sizes are 
equal to the fitting uncertainties. They lie on the SM lines: 
b*o = foo + 0.215ReBo = -0.0280±0.0010. The {6o, ReBo} pa- 
rameters obtained from global fits of KLTK90 [2], BFG97 [1], 
F02a [13] and one of F02b [14] are shown by large open circles 
with both vertical and horizontal error bars. Also shown are 
{6o, ReBo} sets from FG98 [15] (closed squares) and F02a [13] 
(small open circles) , without conversion from bo to fcg in their 
convention. 
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